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ABSTRACT 


Tolerance limits are formulated for the Maxwell distribution, 
and a table of tolerance limit factors for the upper tolerance 
bound is provided as a function of P (percent of the population 
below the bound), y (confidence level), and n (sample size). 
Values of P, y, and n considered are P = .50, .75, .90, «95, .99; 
Y= .75, .90, .95, .993; n = 2(1)25(5)100(10)200 (50)3C00(100)1000,>. 


While the formulation is sufficiently general to be of use to 
anyone who deals with Maxwell data, examples are restricted to the | 
area of weapon systems analysis. It is in this area that the 
analyst is often confronted with the problem of estimating the 
radius of a sphere which will include 100P% of the future burst 
points from an air burst weapon. Under the assumption that the 
distribution of burst points about the target center is trivariate : 
normal with common standard deviation o in all three directions, 
this development will enable him to attach a confidence statemert 
to che percent of the population encompassed. In particular, it 
will enable him to determine, on the basis of test firings, the 
radius of a sphere which will encompass at least 100P% of the 
future burst points with 100Y% confidence. 
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A measure of dispersion often applied to air burst weapons is the 
SEP (Spherical Probable Error). Tnis parameter is the radius of a 
mtan-centered sphere which includes 50% of the trivariate probability, 
er in terms of a particular weapon, it is the radius of a sphere 
centered on the target center within which 50% of the rounds will burst. 
Its bivariate counterpart for ground burst weapons is, of course, the 
CEP which is the radius of a mean-centered circle which includes 50% 
of the bivariate probability. Estimation of the SEP involves firing a 
sample of n independent rounds at a target arbitrarily placed at the 
center of the Cartesian coordinate system. The radial burst distances 


of these rounds from the target center, denoted by las are then 


recorded and used to compute un estimate of SEP, say SEP, which is 
taken as the radius of a sphere within which 50% cf the future rounds 
from this weapon will burst. 


SBP is only a point estimate of SEP and, as such, it will vary 
from sample to sample. To provide a measure of precision concerning 
the percent of the population encompassed, the probability that a 
sphere of radius SEP will include at least 50% of the future rounds 
is considered, Since this probability is found to be quite low, an 
alternative procedure is suggested through the development of tolerance 
limits for the Maxwell distribution, This alternative procedure enables 
one to increase this probability (or confidence) to more reasonable 
levels not only for 50% of the trivariate probability but also for 75%, 
90%, 95%, and 99%. The development is an extension of Thomas et al. 
(1973) for the CEP and is sufficiently general to be useful to those 
other than the weapon systems analyst. 


DERIVATION OF THE MAXWELL DISTRIBUTION 


The use of the SEP as a measure of dispersion requires the assunp- 
tion that the burst points of a weapon about the target center are dis- 
tributed according to the uncorrelated trivariate normal distribution 


with common variance o in all three directions (sometimes referred to 
as the spherical normal distribution). If one lets the trivariate 
random variable (X,Y,Z) designate the miss distance of a burst point 
from the target center (arbitrarily placed at the origin) in the x, y, 
and z directions respectively, thie density is given by 


-3/2 = (x? 4y 7427) /207 
£(x,y,z) = (2n02) e , MO ¥y,ZSS, 
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Consider now the distribution of the radial burst distance from the 


target center, i.e., the distribution of R = (x? + y? + 27%, _cang- 


forming to spherical coordinates by letting 
X= R sin 9 cos @ 
Y = R sin 9g sin 8 
Z=Reos 0, 


it is easy to show that the density of the trivariate random variable 
(2,9,°) is 


Fin 
73/2 .-r°/2o° 12 (2) 


g(r,8,9) = (2n¢”) r sin @ 


for r >0, 6< 8 < 2m, and O< $<. The marginal density of R is 
now obtained by integrating g(r,$,°) over the entire range of both 8 
and @. It turns out that the density of R is 


yo 
h(r) = (2/ny2 (27/0) et fe F r>0o. (3) 


Tris density is comnonly referred to as the Maxwell density and repre- 
sents the distribution of the radial burst distances from the target 
center under the spherical normality assumption. (See, for example, 
Lindgren (1968).) Its counterpart for ground burst weapons is referred 
to as the Rayleigh distribution and is the distribution of the radial 
miss distances in the ground plane. 


Using the density in (3), one can establish the relationship between 
SEP and ao by solving 


SEP 
P(R< SEP) =f h(r)dr = .5 (4) 
0 


for SEP in terms of o. Since h(r) is not integrable in clos2d fom, 
numerical procedures can be employed to show that SEP = 1.52820, Hence, 
if the population standard deviation go were known for a given weapon-to- 
target range, a sphere of radius 1.53820 centered on the target center 
would encompass 50% of the trivariate probability or 50% of the future 
burst points from thisteapon under similar conditions. (In the bivariate 
analogy fer ground burst weapons, it is well known that “EP = 1.1774c,) 


a 


Unfortunately, o and hence SEP are never known and imust be cetimated 
from test firings. To esclmate o, n rounds are fired at a target at a 
given range and the distances of the burst points from the target center 
are recorded. If these distances are designated as Tyo nt is easy to 
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show that the maximum likelihacd estimate for c is given by 


th 


- Di ae 
=a { 
oe lg r,/3n] (5) 
and thus the maximum likelihood estimate for SEP is given by 


SEP = 1.53826. (6) 


Suppose now that cone estimates SEP with SEP as given in (6) above. 

One cannot state that a sphere oi radius SEP will encompass 50% of the 
future pounds from this weapon (under similar conditions) with certainty 
since SEP is a continuous random variable which varies from sample to 
sample, However, as a measure of precision, one could consider the 
probability that a sphere of radius SEP encompasses at least 50% of 

the trivariate probability (at least 50% of the future burst points). 
This probability will be explored in the next secticn. 


PROBABILISTIC DEVELOPMENT 


Thais section concerns the development of an expression for the 


probability that a sphere 2f radius SEP = 1.53826 encompassing at least 
50% of the trivariate probability However, to keep the development 
general, this will be accomplish. by obtaining an expression for the 


probability that a sphere of radius kG encompasses at least 1002% of 
the trivariate probatility. The mathematical expression fox this pro- 
bability is 


ko tia 2 
prob {f  (2/m)® (7/08) a [2° ar 2 P} = y (7) 
ce] 


whcre Y is the quantity sought, Letting y = =*/o" in the integral in 
(7), one obtains 


«6? Io* -4 & a /2 
Prob {J (27) * y*e aus dy 2P}=y. (8) 
° 


The integrand in (8) is tecognized ar a chi-square density with three 
degrees of freedom, Hence (8) car ve written as 


“232 
Prob {F(—2-) 2 P} = Y (9) 
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where F denotes the cumulative distribution function for the chi-square 
with three cegrees of frecdom. Since F is a one-to-one function, (9) 
can be expressed by 


Fl o@j}ey (19) 


Prob ie 5 
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or 


ft 


2-1 
Prob (67 2 oF); 2yY. (11) 
" 


Note thet F_ 1 (py is the 100P percentage point of 5 chi-squeve density 
with three degrees of freedom. Tabular values of pol (P) are available 
from -srious sources, mest notably from Biometrika Tables for Statisti- 
ci-ns, Vol, I where they are presented to six significant digits for a 
wile range of values of P. For P = .50, the value of P under present 


consideration, F590) = 2.36597. 
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Using the maximum likelihood estimator for 6 as given in (5), it 
is easily shown that the density of W = ° is given by 
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3n Ww e 
fe) =, wrod. (12) 
P(3n/2) 22772 g3 
Tnerefore, equation (11) can be expressed by 
v 
J f@w) dw=l-y (13) 
: . 
where £(w) is given in (12) and v = oF “hay in?, It appears from (13) 
that y, the probability urder question, is a function of the unknown 
perameter go. A simple tre eee reveals that it is not. Letting 
zs w/o* so that w = o2z and dw= 242, equation (13) becomes 
Pi 3n 3n/2 _3n/2 =1 _-3nz 
eal ee ge dz=il-y (4) 
d 
° P(3n/2) 2°n/? 
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where v’ = Fl py se? Equation (14) is free of the unknown parameter ¢ 
so that the probability y is a function of only P (% of the population), 
k (multiplying constant for G), and n (sample size). 
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Recall now the question posed at the end of the last section, 
namely, with what probability (or confidence) can one state that a 


sphere of radius SEP = 1.53826 will encompass at least 50% of the 
trivariate probability. This probability can be obtained from 
equation (14) by setting k = 1.5382, P = .50, and solving for y 

for various values of n. This equation was scived with k and P set 

as above for sample sizes of n = 2(1)25(5)100(10)200 (50)306(100)160L,>. 
The results are set out in Tables 1 and reveal that this probability or 
confidence is quite low, i.e., leas than .50 unless one has an infinite 
sample size (tantamount to having complete knowledge about the unknown 
paramete- o). To increase this confidence to more reasonable levels, 
it is clear that one must increase the multiplying constant above 
1.5382, This will be discussed next through the development of 
tolerance limits for the Maxwell distribution. 


UPPER TOLERANCE BOUND FOR THE MAXWELL DISTRIBUTION 


Making confidence or probability statemeats concerning the percent 
of the population which lies belew a estimate of tne SEP involves the 
concept of an upper tolerance bound. In the more general sense, an 
upper tolerance bound, U{P,y), is a point defined such that at least 
100P% cf the population i*es below it with LOOYZ confidence. (See, 
for example, Bowker and Lieberman (1972), Proschan (1953), or Thomas, 
et al. (1973).). U(P,y) is constrected as a function of the estimate(s) 
of unknown population parameter(s) based on a random sample from the 
population in question. For the case at hand in wiich one is sampling 
from a trivariate normal distribution with cocrson variance, there is 
only one unknown parameter, i.c., o, and the upper tolerance bound 
will be formulated as a function of the estimate of this parameter. 
(Since oniy the radial distance between burst point and target is under 
consideration, this can also be viewed as saepling from the Maxwell 
distribution with parameter o. See equation (3).) 


The maximum likelihood estimator for ¢ as given in (5) can be 
shown to be a sufficient estimator for c, so that the upper tolerance 
bound should be of the form U(P,y) = k(P,y,n)o. The constant k(P,v,n) 
(tolerance iimit factor) is to be determined such that one is 16 
fident that at least 100PZ% of the population lies below U(P,y). 
value of P = .50 is of primary concern here, other values could 
of interest so the general notation will be used. Deleting the argu- 
ments for notational simplicity, &(P,v,n) ia scught stch that 
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te aw 2 4» 
or 2 - “ é 
Prob {J Qn)? (r2/0°) e* [20 dr 2P}=y. (15) 
S 
This is, of course, precisely equation (7) of the last section vhich 
reduced to equation (14). While equation (14) was originally derived 
for the purpose of evaluating y for fixed values of P, k, and n, it 


can also be used to evaluate k for fixed values of P, y, and n. The 
value of k so obtained would be an exacc tolerance limit facter which, 
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Table 1 
PROBABILITY THAT AT LEAST 50% OF POPULATION LIES WITHIN SEP 
—=—_e a He OU EULATION LIES WITHIN SEP 
a Y a Y 
2 24232 60 4860 
3 +4373 65 4865 
4 04457 70 4870 
5 24514 75 4875 
6 4556 80 4879 
7 24589 85 4882 
8 4616 90 «4886 
9 4638 95 24889 
10 4656 100 4891 
1l 4672 110 ,4896 
12 4686 120 4901 
13 4699 130 -4905 
14 4710 140 4908 
15 4720 150 4911 
16 4728 160 4914 
17 04737 170 4917 
18 24744 180 4919 
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21 4763 250 24931 
22 4768 300 «4937 
23 4774 400 +4946 
24 04776 500 «4951 
25 «4783 600 «4956 
30 «4802 700 04959 
35 «4816 809 ~4962 
40 4828 900 +4964 
45 +4838 10C0 »4966 
50 +4846 % 1.0000 
35 04854 
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when multiplied times 6, would provide the upper’ tolerance bound 
U(P,y). As aforementioned, the interpretation of U(P,y) is that one 
is 100y% confident that at least 100P% of the population lies within 
a evhere of radius U(P,Y). For a particular weapon at a specified 
m-to-target range, this means that one is 100y% confident that 

least 100P% of the future rounds from this weapon, fired under 
simiiar conditions, will burst within a sphere of radius U(P,¥) 
centered on the target center. 


To utilize the above concept, it is necessary for one to have a 
table of tolerance limit factors at his disposal for reasonable 
values of P, y, and n. Hence, equation (14) was solved for k for 
P= 50, .75, .90, .95, .99; y = .75, .90, .95, .99; and n = 2(1) 

25 (5)100(10)200(50}300(100)1000,%. The solutions were obtained 
using Simpson's integration rule and successive binary cuts begin- 
ning with an appropriate starting value for the upper limit v’, 

The computations were performed on the CDC 6700 at the Naval Weapons 
Laboratory. Tolerances were set to provide an accuracy in k of four 


decimal digits; the values of k are set out in Table 2 for the above 
listed values of P, Y, and n. 


As an example of employing this procedure, suppose eight rounds 
are fired at a target to obtain the radius of a sphere about the target 
center within which at least 50% of the future rounds (under similar 
conditions) from this weapon will burst with 95% confidence. The 
radial miss distances from the target center are shown below. All 
measurements are in feet. 
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201.92 
52.98 
210.10 
120.43 
48.17 
96.40 
85.84 
104.71 
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It can be verified that itt “° = 132,169.9603 and hence that 


i $m 


o = 74,21 


and 


SEP = 1.53820 = 114.15 P 
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From Table 1, it is seen that one is unly 46% confident that a sphere 


of radius SEP = 114,15 will include at least 50% of the future bursts 
from this weapon. To increase the confidence, or probability to 95% as 
specified, one refers to Table 2 under P = .50, y = .95, and n = 8 to 
find the tolerance limit factor k (.50, .95, 8) = 2.0250 (vice 1,5382). 


This is then multiplied times G to cbtain U(.50, .95) = k(.50, .95, 8) 


& * (2.0250)(74.21) = 150.28. Hence, a sphere of radius 150.28 feet 
about the target center will include at least 50% of the future bursts 
fron thic weapon under similar conditions with 95% confidence. Should 
the experimenter want to increase the confidence to 99%, he would refer 
to Table 2 under P = .50, y = .99, n = 8 to find k(.50, .99, 8) = 
2.2870. This, when multiplied times 6, provides U(.50, .99) = 169,72 
feet. Finait,, should one want the radius of a sphere which would 
include 95% of the bursts from this weapon (vice 50%) under similar con- 
ditions with 95% confidence, he simpiy refers to Table 2 to find k(.95, 
.95, 8} = 3.6802. This is then multiplied times @ to obtain U(. 95, .95) 
= 273.11 feet. 


GONCLUS IONS 


The tolerance limit factor, k(P,v,n), derived in this report can 
also be expressed as a function of chi-square oat points. Refer- 


ring to equation (11), it can be shown that 3nd" Io" ig distributed 
according to the chi-square distribution with 3n degrees of freedon. 
Hence, one could write equation (11) as 


Prob ea = aor = ¥ 


where, as before, the arguments of k have been deleted. In order for 
equation (16) to be satisfied, one must have 3nF™) (Pp) /? = Su ey or 
3 


wt 
3nF ~(P 
BaF, 


%n,1-y 


k= 


where Ga, 1-¥ is the 100 G-y) percentage point of a ee density 
with 3n degrees of freedom. Utilizing thie notation, F 1 (py can be ex- 
pressed es 5 P Bo that k(P,y,n) can be written as 


(16) 


(17) 


(18) 


Figen is ain ht biwaseamvennrewemine ome = 


‘Hg -aoaogenet” ' 


Hence, should one need a value 


of k not tabulated 1 


be obtained by refarring to a table of chi-square p 


computing k as in (18) above, 


ri 


n this report, it can 
ercentage points and 
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